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The spin orientation of a magnetic dopant in a zincblende semiconductor strongly influences the 
spatial structure of an acceptor state bound to the dopant. The acceptor state has a roughly oblate 
shape with the short axis aligned with the dopant's core spin. For a Mn dopant in GaAs the local 
density of states at a site 8 A away from the dopant can change by as much by 90% when the Mn 
spin orientation changes. These changes in the local density of states could be probed by scanning 
tunneling microscopy to infer the magnetic dopant's spin orientation. 
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Semiconductors doped with magnetic atoms have at- 
tracted much interest in recent years because of their 
potential applications in spintronic and quantum infor- 
mation technology. The magnetic dopants introduce, 
in addition to magnetic moments, spin-polarized accep- 
tor states into the semiconductor hosti^*^*^*^ Recent ad- 
vances in scanning tunneling microscopy (STM) on GaAs 
have led to a better understanding of the anisotropic 
shape of these acceptor states. ^'^'^"'^^ The principal de- 
terminant of the anisotropic shape is the cubic symmetry 
of the lattice, and not the spin-orbit interaction, as the 
characteristic anisotropic shape is predicted to form even 
in semiconductors with negligible spin-orbit interaction.^ 
However, the spin-orbit interaction does partly corre- 
late the degree of acceptor state anisotropy with the 
spin orientation of the magnetic dopant This corre- 
lation suggests the possibility of detecting the spin ori- 
entation of a magnetic dopant with a purely nonmag- 
netic probe. A similar phenomenon has been found in 
metallic magnetic systems. The spin-orbit interaction 
in iron films mixes the d bands to yield a few percent 
difference in the LDOS for nanoscale magnetic domains 
(out-of-plane magnetization) and domain walls (in-plane 
magnetization) A coupling between spin and orbital 
degrees of freedom in quantum dots has been predicted 
to generate a spin-dependent electric field, although the 
size of that effect has been estimated to be very smallii^ A 
spin-orientation-dependent LDOS has yet to be observed 
in semiconductor systems. 

Here we describe calculations of the dependence of the 
LDOS near a Mn dopant in GaAs on the Mn spin orien- 
tation. The orientation-dependence of the LDOS ranges 
from very high 90%) for tunneling into the acceptor 
state, to quite small (~< 5%) for tunneling into contin- 
uum valence states far from the band edge. Tunneling 
into the continuum valence states resembles the qualita- 
tive behavior seen in iron films;^^ the LDOS orientation- 
dependence is largest near critical points of the electronic 
structure: near the valence band edge for GaAs:Mn, and 
near an avoided level crossing for iron. Tunneling into 
discrete states such as the acceptor level, however, can 
be highly selective of spin and orbital character, as local- 
ized states with different orbital characters have narrow 




FIG. 1: (a) Atomic structure near a substitutional Mn dopant 
(blue) in the GaAs lattice (red atoms are As). The As atoms 
are labeled by SI, S2, S3, and S4. (b-d) Contour surfaces of 
the LDOS of the acceptor level at 10% of the peak value at 
the Mn site. The Mn spin is aligned with the (b) [001], (c) 
[llO] or (d) [111] axis of the GaAs lattice. The symmetry is 
(b) D2d, (c) C2v or (d) Csv The LDOS at each atomic site 
is spatially distributed according to a normalized Gaussian 
with a 2.5A width. The box outlines are aligned with the 
cubic lattice and have widths in units of the lattice constant 
(a = 5.65A). 



enough linewidths not to be spectrally distinguishable. 
Electron spin resonance (ESR) of the Mn spins has been 
observec^ in GaMnAs; thus we predict that the ESR of 
a single Mn spin could be detected using a nonmagnetic 
scanning tunneling microscope. 

The results presented in this paper are calculated from 
a sp^ tight-binding model, including on-site spin-orbit 
coupling, for an isolated Mn acceptor state in GaAs>i^ 
The Mn spin is assumed to be a classical spin aligned 
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in a certain direction, for example by applying an ex- 
ternal magnetic field. The Mn spin degree-of-freedom is 
treated as an effective spin-dependent potential acting on 
the valence electrons of GaAs. Motivated by a p-d hy- 
bridization model,^ our Mn potential has spin-dependent 
matrix elements at the four nearest-neighbor sites. The 
Hamiltonian with a single Mn dopant takes the form 



H = 



(1) 



where Hq is the sp^ tight-binding Hamiltonian for bulk 
GaAs with only nearest-neighbor hopping,- ^ and Vg is 
the effective potential due to the Mn dopant with its 
spin pointing in the direction S. When the Mn spin is 
aligned with the [001] crystal axis, the potential is 



V[001] = Vn >S> 1 + Vm IS> a3 



(2) 



where Vn and Vm are the spatial part (including both 
lattice sites and atomic orbitals) of the nonmagnetic and 
magnetic potential, and tTj's are Pauli spin matrices. In 
our model, Vn is non-zero at the Mn and the four nearest- 
neighbor As sites, and Vm is only non-zero at the four 
nearest-neighbor As sites. The non-zero potential ma- 
trix elements at the neighboring As sites come from the 
hybridization among the Mn d orbitals and the sp^ hy- 
brids of GaAs. The potential V^ for the Mn spin pointing 
in a general direction can then be obtained through a ro- 
tation in the spin space. 
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6 is the angle between the vector S and the [001] axis, and 
n is a unit vector pointing to the direction of [001] x S. 
We find that adding the Coulomb potential out to the 5th 
nearest neighbor to the impurity (and correspondingly 
reducing Vm to keep the acceptor binding energy fixed) 
does not significantly alter the anisotropy of the acceptor 
state we describe in this paper and that is visible in Figs. 
1-4. 

Within this sp'^ model, the energy spectrum of the sys- 
tem is independent of the Mn spin orientation. This 
is a direct consequence of the ability to transform two 
Hamiltonians with different Mn spin orientations into 
each other with a unitary transformation, which wc now 
derive. The basis states (s and p orbitals) are irreducible 
representations of the Td group and of the rotation group 
SO (3). As a result, the Mn potential is invariant under 
a rotation in real space. 



Vs = Ulv^Un 



(5) 



where Uu = U^,^ igj C//?;. The unitary transformations 
and U° correspond to the rotations of the lattice and 
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FIG. 2; LDOS spectra at the sites (a) SI or S2, and at the 
sites (b) S3 or S4 as indicated in Fig. The solid, dotted 
and dashed lines correspond to Mn spin orientations along 
[001], [110] and [110] respectively. Zero energy is the valence 
band maximum. The energy broadening factor for the LDOS 
is 5 meV. (c) The spectral weight at 113 meV as a function 
of Mn spin orientation for the four As sites. 



of the atomic orbitals respectively. Note that the homo- 
geneous Hamiltonian iJg including the spin-orbit inter- 
action is invariant under a full rotation. 



where 



7?o = U^H^U 



(6) 



(7) 



With Eqs. © and the rotation of the Mn spin can 
be achieved via rotating the whole lattice in the opposite 
direction. 



Ho + V. = U^{Ho + Vooi])U 



(8) 



if = [/!. 

Although the energy spectrum (and thus the DOS) of 
the magnetic dopant in a zincblende or diamond symme- 
try semiconductor is independent of dopant spin orienta- 
tion, the LDOS can vary considerably. The symmetry of 
the LDOS would be tetrahedral (T^) if there were no spin- 
orbit interaction and also no Jahn- Teller distortion. For 
a Mn dopant substituting for a Ga atom no Jahn-Teller 
distortion is seen.^ The actual symmetry of the LDOS 
for GaAs:Mn, however, is lowered by the spin-orbit in- 
teraction, and the resulting symmetry depends on the 
Mn spin orientation as shown in Fig. ^ The LDOS is 
calculated as the imaginary part of the Green's function, 
LODS(R) = -TT-iJm [trG^(R,R)]. 

To understand how the spatial symmetry and spin ori- 
entation are correlated we consider the angular momen- 
tum character of the acceptor states. The acceptor states 
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found in our model are spin-polarized antiparallel to the 
Mn core spin due to the exchange coupling between 
the GaAs valence states and Mn d statesi^ Only one ac- 
ceptor state can be occupied by a hole, as the Coulomb 
interaction prevents two holes from binding to the Mn 
dopant. The degeneracy of the singly occupied accep- 
tor states with different orbital angular momentum is 
lifted by the spin-orbit interaction, and the lowest energy 
configuration has an orbital angular momentum aligned 
antiparallel to the Mn core spin, leading to a composite 
J = 1 spin associated with the Mn dopant. As an 7^ 
state of an atom centered at r = has a nodal line along 
X = y = 0, the probability density of the acceptor state 
should be small along the line extending from the dopant 
parallel to the core spin direction. Thus a contour sur- 
face of the LDOS at the lowest acceptor level (shown in 
Fig. ^ has an approximately oblate shape with the short 
axis aligned with the Mn core spin. The shape of the level 
does not depend on which way the Mn core spin is point- 
ing along the axis — only on the orientation of the axis 
itself. One consequence of this shape is anisotropic spin- 
spin interaction;-^^i^^i^''i^^i^^i^°i^^i^^ in which the overlap 
is larger when the two Mn spins are perpendicular to 
the axis that joins them. Figure shows three examples 
for which the Mn spin is oriented along one of the high 
symmetry axes of a cubic crystal. When the Mn spin 
is aligned with the (b) [001], (c) [110] or (d) [111] axis, 
the symmetry of the state is lowered from Td to (b) D2d, 
(c) C2V and (d) C^v . Although the spin direction is indi- 
cated with an arrow, the panels in Fig. would look the 
same if the spin direction is reversed, so for example the 
acceptor state would appear as Fig. ^b) for the Mn spin 
pointing in the [001] or [OOT] direction. 

Figure El shows the energy spectrum near the valence 
band edge and the LDOS spectra at the four nearest- 
neighbor As sites S1-S4. These sites are indicated by 
their labels in Fig. QJa). For a Mn at the origin, SI is 
at (f , -f , f ), S2 at (-f , f , f ), S3 at (f , f , -f ) and S4 
at (— f,— fj— f), where a is the lattice constant. As 
pointed out above, for all orientations of the Mn spin 
the energies of the peaks in the spectra are the same. 
The Mn dopant itself sits at a high symmetry location, 
and thus the LDOS at the Mn site is independent of 
the spin orientation. The spectral weight of the acceptor 
state at the four nearest-neighbor As sites evolves as the 
Mn spin rotates. Results are shown in Fig. ^c) for SI 
(solid hue), S2 (dashed line), S3 (dot-dashed line) and S4 
(double-dot-dashed Hue). From [100] to [001] the spectra 
at S3 (S4) is the same as S2 (SI). From [001] to [110] 
the spectra at S2 is the same as SI. The orbital angu- 
lar momentum aligns with the Mn spin, and thus the 
spectral weight tends to increase when the site is away 
from the spatial line drawn through the Mn dopant and 
parallel to the Mn spin. Conversely the spectral weight 
tends to decrease when the site is close to this spatial 
line. The variation among the nearest-neighbor As sites 
is greater than a factor of two. For sites farther from 
the Mn dopant, the dependence of the LDOS on Mn spin 



orientation can be much greater. The variability of the 
LDOS of the continuum valence band states is consider- 
ably less; the LDOS for all three spin orientations shown 
in Fig. 121 at the valence maximum varies by ^ 40%, but 
is the same within a percent below the split-off energy 
(-350 meV). 

To show the LDOS variability for different Mn spin 
orientations as a function of distance, we plot in Fig. [31 
the probability density of the Mn acceptor state at the 
Ga sites along the [110] direction, which is usually the 
surface normal of a cleaved GaAs sample. The back- 
ground LDOS has been subtracted off so that the be- 
havior at large distances is visible, and the probability 
density has been normalized to unity at the Mn site. 
The background LDOS for the 5 meV energy linewidth 
is about 10"'^ eV""'^. A clear feature of Fig.Ois that the 
probability density drops significantly when the Mn spin 
is oriented parallel to [110] rather than perpendicular to 
[110]. Once the distance is larger than about 4 atomic 
layers (8 A from the Mn site), the probability density 
changes by roughly an order of magnitude as the Mn 
spin orientation is changed from parallel to perpendicu- 
lar to [110]. At large distances from the Mn site we ex- 
pect the wave function follows a simple exponential tail 
with a decay length we estimate based on the effective 
masses of the heavy and light holes and on the accep- 
tor level energy. In our tight-binding model, the effective 
masses along [110] for heavy and light holes are about 
"^hh = 0.5mo and mih = 0.12too. For the acceptor level, 
the binding energy i? = 113 meV and the decay length 
A = n/{2m^gE)'^^/^^ ~ 13 A, where the averaged effec- 
tive mass is rricg = 2mhh'7iih/("^hh + "m-ih)- The visible 
LDOS variability will be significantly reduced at distant 
sites due to the presence of a background LDOS from the 
valence band tail. 

Figure^lshows the field of view at the (110) plane four 
atomic layers from the Mn layer. This is the image most 
relevant to detecting the ESR of a single Mn dopant in 
GaAs. The procedure would be to cleave the GaAs along 
the (110) plane, and locate an isolated Mn a few layers 
down from the surface (here we have chosen four layers, 
but as shown in Fig. any layer would do except the sur- 
face layer). At the central Ga site on the surface plane, 
here located SA from the Mn site in the material below, 
there is approximately a 90% change in LDOS when the 
Mn spin switches from parallel to the surface normal to 
perpendicular to the surface normal. The LDOS also 
differs for the two Mn orientations perpendicular to the 
surface normal, although that difference is considerably 
less (15%). The overall shape of the anisotropic acceptor 
state is very similar for all three spin orientations chosen 
for Fig. 01a-c) , as the gross shape is governed by the cu- 
bic symmetry of the crystal,^ but the degree of anisotropy 
(appearing as the magnitude of the state's LDOS along 
this spatial slice) varies by up to 90%. These characteris- 
tics suggest that atomic-scale resolution is not necessary 
to distinguish the Mn spin axis from measurements such 
as those suggested by Fig. 0] 
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orbit interaction. Calculations of the Mn acceptor state 
were performed within a tight-binding model including 
spin-orbit interaction and p-d hybridization. The accci: 
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FIG. 3: The probability density of the Mn acceptor state 
at the Ga sites along the [110] direction {z is the distance 
from the Mn dopant in units of a). The probability density 
drops significantly when the Mn spin is oriented in the [110] 
direction (square) rather than in the perpendicular directions, 
[001] (triangle) and [110] (circle). The probability density is 
normalized to unity at the Mn site. The dashed lines show 

„-2z/A 



The differences in LDOS shown in Fig. 01 suggest that 
a repetitive preparation of the Mn dopant in a particular 
spin state (using ESR techniques) and measurement of 
the spin orientation at a time delay (using the nonmag- 
netic STM measurement) would be able to produce an 
oscillatory LDOS. The temporal frequency of the LDOS 
created by the Mn spin's precession would be determined 
by the magnetic field strength and g factor of the Mn 
dopant (measured^ to be 2.77), and the spatial struc- 
ture could be determined from Fig. 0] Based on the 
linewidth in bulk ESR measurements done at 9.4 GHa^ 
(< 0.05 Tesla), the spin coherence time of a single Mn 
spin is estimated to be longer than 0.5 ns. Ultrafast 
STM detection at 50 GHz has been achievediS If the 
static magnetic field of an ESR apparatus were directed 
parallel to [001], then precession would involve an oscil- 
lation between Fig. ^fb) and Fig. QIc). As the Mn spin 
has J = 1, pulses of oscillating magnetic fields could be 
used to manipulate the populations of the Jz = 1, 0, —1 
eigenstates. For J[ooi] = ±1 the LDOS would appear as 
Fig.^a). The orbital wavefunctions transform according 
to T2 in the tetrahedral crystal, so the spatial structure of 
|V'J.=oP can be written as |i/'j,=i|^ [V^j^^iP - IV'J.=iP- 
Thus the spatial structure of the J[ooi] = state would 
be Fig. Efa) subtracted from the sum of Fig. EJb) and 
Fig. ^c). The difference between this J[ooi] ~ state 



and the J[ooi] = ±1 states is shown in Fig. gjd). 

We have shown how the anisotropic spatial shape of the 
Mn acceptor state, which principally originates from the 
cubic symmetry of the GaAs lattice, is influenced by spin- 




FIG. 4: Cross-sectional view of the LDOS on the (110) plane 
when the Mn dopant is four atomic layers (corresponding to 
2 = 2''^^^^ a in Fig.EJ below the viewing layer. The Mn spin is 
orientated along (a) [001] (b) [lIO] or (c) [110]. The absolute 
difi'erence between (b)-|-(c) — (a) and (a) is shown in (d). The 
shape of the apparent feature is largely the same for (a-c), 
but the amplitude differs considerably: the spectral weight 
at the center changes by ~ 90% between (b) and (c), and by 
^ 15% between (a) and (b). The LDOS at each atomic site is 
spatially distributed according to a normalized Gaussian with 
a 2A width. 



tor state has an approximately oblate shape with the 
short axis aligned with the Mn core spin. The difference 
in the LDOS for Mn spins oriented in two different direc- 
tions at a site located SA away from the Mn dopant can 
be as high as 90%. We suggest that the high visibility 
of the spin orientation in the LDOS for Mn spins can be 
probed by scanning tunneling microscopy, and this can 
lead to a nonmagnetic detection scheme for single-spin 
ESR. 
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